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The method of calculating the free energy and thermodynamic characteristics of the 
classical n-vector three-dimensional (3D) magnetic model at the microscopic level without 
any adjustable parameters is proposed. Mathematical description is perfomed using the 
collective variables (CV) method in the framework of the p 4 model approximation. The 
exponentially decreasing function of the distance between the particles situated at the N 
sites of a simple cubic lattice is used as the interaction potential. Explicit and rigorous 
analytical expressions for entropy,internal energy, specific heat near the phase transition 
point as functions of the temperature are obtained. The dependence of the amplitudes of 
the thermodynamic characteristics of the system for T > T c and T < T c on the microscopic 
parameters of the interaction potential are studied for the cases n — 1,2,3 and n — > oo. 
The obtained results provide the basis for accurate analysis of the critical behaviour in 
three dimensions including the nonuniversal characteristics of the system. 



I. INTRODUCTION. 

Investigating the behaviour of the real three-dimensional systems near the phase transition (PT) point 
is one of the most important problems of condensed matter physics. In the present paper we propose a 
theoretical description scheme of the critical behaviour of the classical n-vector magnetic model in three 
dimensions on the microscopic level without any adjustable parameters. The description is based on the 
original method of calculating the thermodynamic and structural characteristics of 3D model systems 
near the PT point, which is known as the collective variables (CV) method Q|, 0. Within the frame of 
an uniform scheme this method allowed to obtain universal characteristics of the system such as critical 
exponents and critical amplitude ratios for thermodynamic functions above and below the critical point as 
well as non-universal ones. Starting from the initial principles we calculate the critical temperature and 
obtain explicit analytical expressions for basic thermodynamic characteristics of the system as functions 
of the temperature and the n component number of the model. For the first time the dependence of 
the thermodynamic characteristics of the system on the microscopic parameters of the initial interaction 
potential is investigated. The investigation of the non- universal properties of the system near the PT point 
require the use of non-Gaussian distributions of the spin density fluctuations. The suggested calculation of 
the partition function is based on the use of the simplest non-Gaussian measure density, i.t. the p 4 model 

@ - n 

We have devoted this paper to the widely studied classical n-vector magnetic model m on three- 
dimensional simple cubic lattice, which is also known as the Heisenberg classical 0(n) spin model or, in 
field theoretic language, as the lattice 0(n) nonlinear a model. The investigation of the critical behaviour 
of the classical n-vector model and its partial cases was conducted in the frame of different methods such as 
high- and low- temperature series, the field theory, semi-microscopic scaling-field theory and Monte-Carlo 
calculations. The main attention in these works was given to the investigation of the universal properties 
of the system such as critical exponents ||- jl2| and some combinations of the critical amplitudes of ther- 
modynamic functions Ol- pq]. Besides, the equation of state of Ising system was obtained to order e 2 by 
Avdejeva and Migdal P7| , Bresin, Wallace and Wilson p8| . Later these results were generalized for the 
case of the n-vestor model p8| . 

Some important results were also obtained for calculating the thermodynamic functions near the critical 
point. One of the first works in this sphere is a paper by Wegner (the so-called Wegner's expansion) 
p9|, suggesting the expression for the free energy with "irrelevant" operators in Wilson approach taken 



into account. The works by Fisher and Aharony |30|, Nicoll and Albright pN, Nelson |32| were also 



devoted to receiving of crossover scaling functions for temperatures T > T c in zero external field near 
four dimensions. In 1974, Riedel and Wegner |?3| developed a numerical technique, termed the scaling- 
field method, for obtaining crossover scaling functions for the free energy and susceptibility. Crossover 
functions rather than power laws [B4|- [B6J have described the nonasymptotic region between criticality 
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and the noncritical "background". In the frame of the massive field theory by Bagnuls and Bervillier 
jTof , |37|] the nonasymptotic critical behaviour for d = 3 in the disordered phase case was analysed. They 
obtained explicit expressions for the correlation length £, the susceptibility x an d the specific heat C as 
temperature functions in the disordered phase along the critical isohore for one-, two-, three - component 
models. The description of nonasymptotical (though still critical) behaviour was obtained as crossover 
between the Wilson-Fisher (near the critical temperature T c ) and mean-field (very far from T c ) behaviours 
with three adjustable parameters used. 

At the present time the actual task of the critical phenomena physics is elaborating the methods giving 
quantity description of the critical behaviour of the system without using any adjustable parameters. This 
has been demonstrated, for example, in the works of Dohm and coworkers p8[- pf| ] in which calculation of 
the temperature dependence of the thermodynamic characteristics of the system was performed without 
using e - expansion in the frame of the minimal subtraction scheme for the n- vector model in three 
dimensions on the basis of high-order perturbation theory and Borel resummation. In these papers the 
amplitude functions of the susceptibility, the correlation length, and the specific heat above and below T c 
up to two- loop order within the 4> 4 model for the cases n = 1,2,3 were calculated. This approach exploits 
simultaneously the experimental information and simple relation between the specific heat above and below 
T c for defining the effective renormalized static coupling of the model in the terms of the measured specific 
heat. The effective renormalized static coupling determined in such a way allows to obtain expressions for 
other thermodynamic characteristics of the model above and below T c without using additional adjustable 
parameters. Besides, in the recent works by Butera and Comi high- and low- temperature expansions 
for the free energy, the susceptibility and the second correlation moment of the classical n-vector model 
on the simple cubic (sc) and the body centered cubic (bec) lattices were extended to order 21 . This 
research only contains temperature dependence of the thermodynamic characteristics and does not give 
the possibility to describe the functional dependence of basic thermodynamic functions on the microscopic 
parameters of the interaction potential and characteristics of the crystal lattice. All this indicates that 
nonuniversal properties of the 3D system near the phase transition point have not been studied sufficiently 
yet. However, the precise role and significance of a lattice structure and interaction potential parameters 
for the approach to asymptotic critical behaviour still seems open to question. We hope that our explicit 
representations may provide useful benchmarks in studying this question. 

The approach to the investigation of the critical properties of the n- vector magnetic model, suggested 
in |42] ] provides the necessary conditions for our complex approach to the study of the universal and non- 
universal phase transition characteristics. The approach suggested in this paper allows us to perform the 
analysis of the dependence of the thermodynamic characteristics of the n-vector 31? magnetic model in 
the vicinity of the phase transition point as functions of temperature and study their dependence on the 
microscopic parameters of the interaction potential and characteristics of the crystal lattice without any 
adjustable parameters being used. These results are interesting from the point of view of comparing the 
theoretical investigations and experimental data. 

In Sec. II of the present paper we perform the calculation of the partition function of the n-vector model 
using non-Gaussian measure density. The explicit analytical expressions for partial partition functions and 
general recursion relations (RR) between coefficients of the " effective Hamiltonian blocks" which arise 
in that case are obtained in the p 4 -model approximation. We define the scope of application of these 
approximate solutions and show that RR as partial solution have a saddle-type fixed point for all n. We 
perform the calculation of the eigenvectors and eigenvalues of the RG transformation matrix and give the 
results of the investigation of the dependence of the PT temperature on the microscopic parameters of the 
interaction potential and the characteristics of the crystal lattice. 

Sec. Ill is devoted to the calculation of the free energy of the n-vector magnetic model for temperatures 
above and below the phase transition point. The main idea of such calculation lies in considering separately 
the contribution from the critical region (CR) where renormalization group symmetry takes place and the 
region of the long-wavelength fluctuations (LWF) of the spin moment density. It shown that in the case of 
the temperatures T < T c in the region of the LWF the fluctuations are described by non-Gaussian distri- 
bution with negative coefficient at square term. The distributions of the spin moment density fluctuations 
after the selection of the ordering free energy is reduced to the Gaussian distribution. The dependence of 
the coefficients of the complete expression of the free energy on the microscopic parameters of the initial 
interaction potential and characteristics of the crystal lattice is investigated for the cases n = 1,2,3. The 
contributions into the expressions for entropy and specific heat from the CR and LWF region are analyzed. 
It has been shown that considering the contribution of the LWF region satisfies the positiveness of the 
specific heat of the n - vector model and system stability. 



2 



In Sec. IV. the explicit expressions for thermodynamic functions of the model as functions of the tem- 
perature are obtained. The dependence of the critical amplitudes of the thermodynamic functions on 
the microscopic parameters of the interaction potential and the characteristics of the crystall lattice is 
investigated for the cases n = 1, 2, 3. It has been shown that in the case n = 3 (the Heisenberg model) the 
specific heat at the T = T c has the finite value. The dependence of the maximum of the specific heat on 
the microscopic parameters of the interaction potential is examined for the case n — 3. The ratio of the 
critical amplitudes of the termodynamic functions at the temperatures T > T c and T < T c is calculated 
in order to compare the obtained results with the results of the other methods. 



II. CALCULATION OF THE PARTITION FUNCTION AND INVESTIGATION OF THE 

RECURSION RELATIONS. 



The critical behaviour of the different physical systems is characterized by their belonging to a specific 
class of universality which is defined by the dimensionality of the system d and the symmetry of the order 
parameter n. The Stanley model Q is selected as the object of our investigation of the critical phenomena. 
This model describes the system of interacting n-component classical spins localized at the N sites of the 
rf-dimensional crystal lattice. The Stanley model is a generalization of the series of different models. In the 
case n = it is reduced to the task of a self-avoiding walk and is used for describing the polymerization 
phenomena. The cases n = 1,2,3 correspond to the Ising model, the XY - model and the Heisenberg 
model respectively. The boundary case n — » oo is equivalent to the Berlin-Kac spherical model [ fl3"[ for 
which the exact result is known. The Stanley model is described by the Hamiltonian 



r,r' 



R-R 



SrS 



R3R, 



(2.1) 



where Sr 



Sji ) is the n - component classical spin with the length m 
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Spins are localized at the N cites of the d dimensional simple cubic lattice with coordinates R. The 
interaction has an exchange character and can be described by the exponentially decreasing function of 
distance beetwen the particles. 



J 



( R - R' ) 



Aq exp 



R-R 



(2.2) 



where Ao,b are the constant. In the CV /?k = \P\^ > ■••> J representation the partition function of 
the model ([□]) is written as [J, Q 



J\p]{dpi 



(2.3) 



Z = J exp i^/3$(fc)p k p_ k 
where $(fc) is the Fourier-transform of the interaction potential (2.2), the element of the phase space is [] 



{dp*) N = X\dp«\{dp^dpl 



k i 



1 The prime means that the product over k is maintained in the upper half-cube of the Brillouin zone and k ^ 0. 
The wave vector k assumes all the values inside the first Brillouin zone 

ki = 2-KTii/N c, i = 1, .., d, 

where m are integers, 

=^<ni<^-, (N" ) d = N, 
and iV is the number of the particles in the periodicity volume V = Nc d , c is a constant of the simple cubic lattice. 
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and J[p] is the transition Jacobian from the spin variables to the CV. The expression for J[p] is given in 
Appendix 1. In the coordinate CV representation 

p(R) = J2 Pk exp(-ikR) (2.4) 

the expression for J[p] is factored 

J[p] = exp«V) J'[0] J] J[p(R)} (2.5) 

R 

where 

J[ P (R)]=^-|L|p(R)| 2 '. (2.6) 

i>i ^ 



The expressions for an are shown in Appendix 2. The partition function (2.3) contains the contributions 
of two different types. The first is energy contributions 

\ ]T /3$(fc)p k p- k . (2.7) 

k<B 

They are diagonal in p^ representation and connected with the interaction potential. The second is entropy 
contributions 



£>J(p(R)). (2.8) 



i=i 



They are diagonal in p(R) coordinate CV representation. There are two possible approaches to calculating 
(2.3). The first one is using pk variables in the calculation of the partition function of the system. In such 
a case the energy contributions are diagonal and the nondiagonality of the entropy contributions leads to 
the approximate method of calculation. One of such approximations is writing (2.3) in the following form 



Z = j exp(-- («2 - /8$(fc))^_k) 



k>B 

x[l + r ? +ir ? 2 + ...](dp k ) w , (2.9) 



where 

\ - a 2 i 



^ (21)1 

l>2 y 1 kt 



V 1 ^^p kl ,...pk 2! <5 kl+ ... + k 2m . (2.10) 



As rule, such a way of calculation assumes that for the value rj it is sufficient to restrict the consideration 
to one term (I = 2) and in addition it is presupposed that || 

(a 2 -P$(k)) 1 »a 2i . (2.11) 



The second way of the calculation of the partition function (2.3) assumes the use of the p(R) CV repre- 
sentation. In that case the entropy contributions are diagonal (2. 5), (2. 6) and approximation relates to the 
energy contribution (2.7), as it is nondiagonal in the p(R) representation. The essence of such fitting leads 
to a certain approximation of $(fc) (see. 0) which allows to diagonalize the expressions in the exponent 
under the integrand for the partition function 



R J 



^ a21 p(R) 2l )dp(R). (2.12) 



We use such approximation for <fr a pr in which & a pr is a constant for every interval k S {Bi,Bi-\) and 
equals the respective mean value $(fc). The first way corresponds to considering the moments of a certain 
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Gaussian distribution. It allows us to calculate only certain classes of the graphs and does not solve the 



problem of the description of the critical behaviour as such. Besides, that condition fl2.1l|) is too strong in 
the vicinity of the phase transition point. Nevertheless, such approach has its advantages at the expense 
of considering the wave vector dependence of <&(fc). The second way makes it impossible for us to consider 
the dependence of $(fc) on the wave vector k on the boundary of the Brillouin zone.0 Nevetheless, such a 
way is not restricted to Gaussian moments and is based on the use of the non-Gaussian measure density 
. In this case we do not need to perform the summation of various types of infinite series of perturbation 



theory, the certain terms of which tend to infinity with T — > T c . Besides, the condition (2.11) becomes 
optional. In our opinion this condition is the basic barrier in the description of the critical behaviour of 
the three dimensional systems. The calculation of the nonuniversal characteristics of the phase transition 
particulaly the PT temperature T c is connected with the choice of the interaction potential. The Fourier 



transform of the interaction potential (2.2) takes the form 



$(0) 

" (TTW' (2 ' 13) 

where $ (0) = 8itA (|) . The value of $(fc) for the wave vectors similar to the boundary of the Brillouin 
half zone (B = — ) is much small than $(0). In this region of the wave vectors a weak dependence of <!>(£;) 
on the wave vector is observed. This allows to accept the following approximation for $(fc) 

m = f *(0)(l-2^),*<*' 

\ $ = const, B < k < B. y ' 

The coordinate B is obtained from the condition of the applicability of the parabolic approximation for 
$(fc) and equals 

B = (bV2) 1 . (2.15) 



Among the set of the CV pk there are those connected with order parameter. In the case of a model 
with the exponentially decreasing interaction potential (2.2) it is the po variable. The investigation of 
the critical behaviour of this model is largely determined by considering the contribution from the po 
variable in calculating the free energy. The mean value po describes the behaviour of the order parameter. 



Nevetheless, as seen from the (2.3),( Al.l), all the CV pk are in terc onnected, and the contribution of the 
variable po alone cannot be separated in the partition function (2.3). The given task can be accomplished 
in case we use the method suggested in Jl} . Its essence lies in sequential integrating of the pic variables 
with k and the investigating of the functional from the po variable. In the case of the n-vector model 
such calculations are made for the first time. The functional representation of the partition function of 
the n - vector magnetic model considering the interaction <f>(fc) at k close to the boundary of the Brillouin 
zone has the form Q 



Z= j'[0]exp(u o A') J exp -~ ^ d(fc)p k P-k 



k<B' 



CI4 

4! N' 



2J Pki---Pk 4 , 5k 1 + ...+k4 



ki,...,k 4 ,fci<B 

Part of the interaction potential is contained in the coefficient 

d(k) = a 2 - /3$(k). 



(dpi 



(2.16) 



2 The method considering the dependence on the wave vector k in the frame of the non-Gaussian measure density 
was proposed in the work 0. It should be noticed that the r\ correlation function critical exponent in the case of 
Ising model was calculated in these works. 
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Besides, the coefficients a-n have as a small correction ($ < $ (0)) that part of the Fourier-transform of the 
potential $(fc) which corresponds to the values of k in the vicinity of the Brillouin half zone (see Appendix 
2). The basic idea of the calculating ( [2.16| ) lies in a sequential exclusion from the consideration of the 
short-wavelength variables /?k which decribe the behaviour of the effective spin blocks whose size is small 
in comparison with the system correlation length as the function of the temperature. After each of such 
step- by-step exclusion the sizes of effective spin blocks are increased by the factor of s (s> 1). The set of 
the CV pk is divided into subsets. Each of these subsets contains the variables p*k with certain values of 
the wave vectors k. For the I subset we have k € (Bj+i-B;), where -B;+i = and B = B . In each of 
the layers of the CV phase space the value <I>(fc) is replaced by the corresponding mean value @. After 
the layer integration the partition function ( 2.16| ) can be represented as a product of the partial partition 
functions Qi of separate layers 



where 



Z = C l Q Q 1 ...Q l I {dp k f l+1 w l+1 (p), 



(2.17) 



(P) = exp 
Jn,l+1) 



-\ £ d^(k)p k p^ 



k<B l + 1 



A\N l+1 



Pki— Pk 4 <5ki + ...+k4 



ki,...,k 4 ,fci<i 



(2.18) 



C{=^2 {Nl - Nl - x)n , Qo = [Q(u)Q(d )] N , 

Qj^QCdOQ^-i)]*'. (2-19) 

The value Inuji+i (p) corresponds to the I - block Hamiltonian which depends on the Ni variables /?k 
connected with the fluctuation of the spin moment density in the blocks. The values which are part of the 
expressions for the partial partition functions are written in the form 

Q (*) = ^ (4$) 4 u (^ x ) exp (t 



' 9 (n,i-l) \ ~ 

Q(P l _ 1 ) = (2n)-*\^s d ^ x 



xtf(^^,W-i)e*p(^)> (2.20) 

Q N ' ( u ) = j [0] exp Cu N'^ . (2.21) 
For coefficients a^' and a.4 ' the recursion relations (RR) are valid 

a (n,l + l) = fl (rU) + d n ^ Mn ^) > 

4"' i+1) = at' l) s- d E n ( Xl ) , (2.22) 

where d n (-B;+i, B{) corresponds to the mean value in the / layer. Specific functions introduced here have 
the following form 

En ( Xl )=s M ^\, M n ( Xl )=N n ( Xl )-l, N n ( Xl ) = ^M- , (2.23) 

<Pn{Xl) XiU n (Xi) 

where the functions <p n {t), U n (t) are combinations of the parabolic cylinder functions U(a,t) 
<p n (t) = (n+ 2)Ul(t) + 2tU n (t) - 2, U n {t) = j^^jy 
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The variables 



>.0 



d n (Bi + i,Bi) , 



Vl = s 2 U n (xi) 



fn (Xl) 



(2.24) 



serve as arguments, where Ni — N s , d is the dimension of space, I is the layer of integration. The 
obtained expression ( [2.17 ) for the partition function enables us to calculate the free energy of the system 



F= -kTj2^Qi 



(2.25) 



when the explicit analytical expression for the partial partition function Qi ( |2.19|) , (2.20) arc familiar 
to us. But, the calculation of the ( 2.25 ) can only be done if the explicit solutions of the RR ( 2.22 ) are 
obtained as functions of the phase layer number I. Similar investigation of the RR for the n - vector model 
was carried out in series of works p2]] , ]45[] , ptfij] , [ ft7| . The main attention in these papers was given to 
calculating the correlation function critical exponent and other universal characteristics of the model for 
the marginal cases s — > 1 , jl5| , |47) and s ^> 1 J46|. In this work we briefly present the results of the 
investigation of the RR difference form (with s > 1) |44[| which allow us to perform the calculation of the 
non-universal system characteristics. We introduce the following definitions 



» 



s M d n (0), u 



(n) 



with 



d n (B l+1 ,B l ) 



,.2/ 



(2.26) 



(2.27) 



where 



q = P$(0)q, q 



(1 



-{d+2)) 



id- 



(1 



(2.28) 



q is the geometric mean value of fc 2 in the interval 1). As a result, the obtained RR (2.22) have the 
following representation 



(n) 
r i+i 



Q + (r^ + q)N n {xi) 



^+1 = S^U^Enixi). 

The equations ( p. 29 ) as a partial solution have a fixed point 

< = -/„/3<f(0),< = ^(/3$(0)) 2 , 

where 

s 2 (N n (x*) - 1) _ 



(2.29) 



(2.30) 



s 2 N n (x*) 
3 



1 *' 

1 - s- 



(x*)^ {N n (x*)-s- 



(2.31) 



The constants /„ and tp n depend on the component number n of the model and the universal value x* 
which is a solution of the nonlinear equation 



„4+( 



(p n (y(x*)) = ip n (x* 



(2.32) 



The values /„ and ip n for different n and s are presented in Table|. 

The presence of the fixed point in RR ( 2.2S ) allows to write them in the linearized form 
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= 5ft 



(n) 
1 

,(n) 



r 



(2.33) 



In calculating the matrix elements of 5ft matrix we restrict our consideration to linear on (xi — x*) approx- 
imation. As a result the following expressions are obtained for the matrix elements (ij] |j 



Ru = s 2 VSfii, R12 = — (po - Hxx*)(u*) 1/2 , 



R21 



„4-d 



7 3u*wi, R22 = s 



4-d 



w 



LO\X 



The eigenvalues E\ and E% of the matrix 5ft are universal values |44 



49 



E%, 



R22 ± 



(i?H — R22) 2 + 4:Rl2 R21 



1/2 



(2.34) 



(2.35) 



As is obvious from Table[| we have a saddle- type fixed point E\ > 1,E2 < 1 for all values n and s. 
According to E§] , E9| the bigger eigenvalue Ei defines the critical exponent v for the correlation length 



Ir 



In E x 



(2.36) 



Since E\ is a universal value, the critical exponents are also universal and do not depend on the mi- 
croscopic characteristics of the system. As we can see from ( 2.34 ) and ( [2.3E| ), they only depend on the 
dimensionality of space d and the n spin component number. The results of calculating the critical expo- 
nents of the model in the framework of this approach was shown in f5(i[| , fsifl . Thus, in the case d = 3,n. = 3 
the following values v = 0.674, a = —0.021,7 = 1.347 were obtained. These values of the critical expo- 
nents correspond to the p 4 model approximation which gives a good qualitative description of the critical 
bahaviour of the n- vector model. As is known ||, good quantitative results for the critical exponents 
can be obtained in the framework of the p e model approximation. For example, the value of the v critical 



3 It should be pointed out that in (2.34) we introduce the following definition 



/'J = /'.. [ 'U - y ) I /'" = 



1/2/ 



* 2 V(y*) 



■ uj\ = uja{bi — qi), 



a 1 =P 1 y*r 1 , ri=S x -y, P m 



bi = Qiy*ri, 



U n (y* 



d m U n ( yi 



<fin(y*) \ dy\ 



d m Vn ( yi ) 



1/2 p(0) 



-1/2 



4 In the marginal case of the larger n — > 00, the cumulants U21 strive for their limit values 

6 



lim U2 = 1, lim 114 = lim 



n + 2 

lim U21 =0, I = 3,4, .. 



0, 



(in normalizing the spin dimensionality m 2 = n). Considering only the leading by order of the value 1/n members 
in the process of integration of the partition functions enables us to employ the method of g-expansion jl2|, | ^5[ , 
p7| . In this case the values 

U21 i ^ d{B u B') \ 



mi 



{2l)\> I 



9 = (ff4,Se,-) 



are small. It allows obtaining the next relations for the correlation length critical exponent (in the case r\ = 0) 
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exponent in the p 6 model approximation increases and practically remains unchanged with the increase 
of the m order of the p 2m model. Obtaining the eigenvectors of the transform matrix !ft from ( 2.33 ) is an 
essential aspect of the investigation of the RR. They can be represented as follows W , |44| 



where 



W 1 = Wu 

Ri -- 

R 



1 

Ri 



R21 



W 2 = W 2 2 



Ei — R22 

R\2 



E\ — Rn 

R12 
E2 — R22 



E 2 — Rn R 2 i 
The conjugate vectors vi and V2 are written in the form 



Vi = Vn 1 



Ri 



Ei — R 2 



v 2 = V 2 2 



Eo — R 



22 



Ri 



(2.37) 

(2.38) 
(2.39) 

(2.40) 



The normalization conditions wivi = 1, W2V2 — 1 give the relations for obtaining the coefficients wu, va(i = 
1, 2). Proceeding from (|2.33| ), ( ^3?| ) the RR (|^|) can be written in the form 



ciE[+c 2 RE l 2 , 
= u* n + ciRiE[ +c 2 E l 2 , 



M * 
r\ = r ri 

,(«) 



(2.41) 



where ci = c'iWh; C2 = c' 2 u>22; c[, c' 2 — const. 
From the initial conditions at I = 

the following expressions for the coefficients 1 , 2 were found 



(2.42) 



C2 



where D = El „ Rl1 = S 2 El = El E r? . For temperatures similar to T c the coefficients ci and c? can be 
represented as follows 



S* = 0, v=2' at d>4 ' 



92 = 



(n + 8)(l 



2 + 



In s 



, at d < 4, 



where <?* are the fixed points. But, starting from the type of fixed point we have certain restrictions on the s value 
(1 < s < 2). In the limit n — > 00 for d < 4 we obtain 1/ = 1, which is in agreement with the Berlin-Kac spherical 
model Analogical expressions for v were obtained in a series of the works E3], EEI, p7| , E3|. 
5 The relations for obtaining the coefficients Wii, are written in the following form 



wnvn 



1 + 



R12R21 
{E1-R22Y 



1 + 



{Ei - Ru) 
R12R21 



li'221'22 



1 ' 



R12R21 
(E2-Rll) 2 



Cgg - ^22) 2 
-R12-R2I 



9 



Cl (T) = c 1T r/3$ (0)), 
c 2 (T) = c 2T \p* (0)] 2 . 



(2.43) 



Having applied these relations we evaluate the critical region of the temperatures r < r* in which the 
solutions of RR (2.41) are valid. For the critical region to exist it is necessary that the "exit" from this 
region at I — > 1 should not exceed the " entrance" . It means that the value r* is equal to the smaller root 
of the two equations 

c 2 RE 2 = CiEi, ciRiEi = c 2 E 2 . 



Take into account (2.43) we obtain for r* 2 

r * = | c 2T i?^ 2 /3"(0) | 
1 citEi 



c 1T E 1 R 1 



The results of calculating t^ 2 in the case s = 4 are shown in Table||. As c\{T) ~ r, the value cyt can be 
represented in the following approximation 



cit = cife + ci k iT + O (r 2 ) 



where 



Clk 



Clfel 



Cll 



Cl2 



D-\ 



(2.44) 



(2.45) 



[&*(o)r 



ith 



en = 1 -/„- J2Vn /2 . 
C12 = -a 



In accordance with (2.43) for c 2 t we consider only the terms proportional to r 2 

C2T = C 2 fe + TC 2 fcl + r 2 c 2k2 + O (r 3 ) , 
where the following definitions are introduced n 



(2.46) 



(2.47) 



C 2 k 



C 2 kl 



C2fc2 



C23 



C22 



C21 



PMo) \p c $(o)] 

C22 2c 2 l 



D-\ 



PMO) [/3 C $(0)] Z 



D~ 



(2.48) 



C21 



tD- 



The calculation of the partial partition functions (2.19),(2.2C) are connected with the employment of 
common RR (|2.22j ). In the vicinity of the fixed point ( 2.301 ) they can be replaced by the approximate 
relations ( 2.41| ) which are exact for T — T c . The q uestion arises under what conditions the relations (2.41) 
can be made use of instead of common RR (2.22). It is obvious that the system of the relations 



It should be pointed out that 



(n,0) (n,0) „* 1/2 

C21 = a\ ,C 2 2=— a) 2 ' R±ifn i 

C23 = Rl¥>n 2 (1 - fn) - <fn, 

R* = Ri (u n y 1/2 ,R* =R{u* n f /2 
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» 



(2.49) 



I ( n ) * | ^ * 



is a condition of the applicabilit y of the approximate relations ( 2.41 ). 

The magnitudes which form ( [2.4S ) are connected with the value of the n spin component number and 



the value of the interaction number I. The main reason for the deviation of the values 



» 



and 



,(«) 



from 



their fixed values is the availability of the terms proportional to c\E\ in the solutions (2.41). For small 
values I the contribution of these terms is small us compared with r* and w* , as c\ ~ r. But in the case 
T 7^ T c there always exists such a value as I — m T so that the contribution will be of the order r* or 
«*. With I > m T the deviation will become considerable and the equations (2.41) cannot be used for 
the description of common RR (2.22). It is necessary to use the first equation in (2.41) for determining 
the value m T . It is connected with the presence of the R\ coefficient near the term C\E\ in the second 
Proceeding from the assumption that R\ is small value as was shown in the case n = 1 



equation of (2.41 



j3| , we see that the deviation rj from r* will be going on faster than the deviation it ; 
it follows that the value m T may be determined from the equation 



(n) 



' m T + l 



= Sr 



from u*. Hence 



(2.50) 



where 6 is a certain constant value (6 < 1). Since later we will compare the results obtained in the case of 
T > T c with those o btain ed for T < T c , we assume that (5=1 with T <T C and S = -1 with T > T c . The 
condition similar to (p.5C|) was employed in pj.In the first approximation on [Eil E\) mT for m T we have 



In 



lnE t 



+ to - 1, 



(2.51) 



where 



TOo = TO C + TOiT. 

For the coefficients to c and mi the following relation take place 



ln(f n S/c lk ) 



, TOl 



Clfcl 



InEi ciklnEi 
The obtained value of m T defines the point of exit of the system from the critical region 

B m ^=B's- m ^. 

The analysis of the RR shows, that in the vicinity of the PT point two different fluctuation processes 
take place. The first one corresponds to the index values I € (0, m T ) and describes the renormalization 
group symmetry that takes place in the vicinity of the PT point. It is so called critical regime (CR). The 
second one corresponds to the index values I > m T and describes the long-wavelength fluctuations (LWF) 
of the spin and is valid both near and far from the PT point. It is characterized by Gaussian distribution 
with dispersion depending on the availability of the CR. It is so called Gaussian regime for T > T c and 
the inverse Gaussian regime for T < T c . 

The calculation of the PT temperature is a significant moment of the investigation of the critical be- 
haviour of the model. As was shown in the works pl|, ||, the critical temperature is the point where 
Gaussian regime is absent and subsequent relations take place 



lim rll\(T c ) 



I — *oo 



lim r 

l — >-oo 



(») 



(T c ) 



lim w[?j_(T c ) = lim u\ n) (T c ) 

l— *oo I— >oo 



const, 
- const. 



(2.52) 



In accodance with ( 2.41 ), this condition is realized only in the case 

ci(T c )=0. 

With ( [2.42| ) taken ino account the explicit equation for critical temperature was obtained 



(2.53) 



Jn,0) 



l-fn 
p* — 1/2 



>,o) 



PMo) 



(2.54) 
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where R* = R^/u^, and a^ 1 ' ^, a^"' -* are functions of the initial interaction potential (2.14) (53]. In this 
paper we present the rezults of calculating the PT temperature with the interaction potential of the (2.14) 
type. As follows from ( p. 14 ), the obtaining of the concrete calculation results are connected with the choice 



of the value 4>, that is the value $(fc) with k G A. The correction, considering the influence of the Fourier 
transform of the potential in the interval k £ [B ,B), makes the results of the calculation more precise. 
Let us choose the next form for $ 

*= <$(fc))+$ oc , (2.55) 

where 

, , „ f^dk^>(k)k 2 

= rB-ui ■ (2 ' 56) 

J B , dkk z 

The results obtained for the temperature in the limit - — > oo must be in agreement with the mean field 
theory results, i.e. 



Take into account this condition and the equation for the temperature T c (2.54), we obtain the equation 
for defining the $00. The solution of this equation is written as 

$oc = _(" + 2)(/n (2 . 57) 

3n(l-V*) 

The results of calculating the PT temperature of the n- vector model in the case m? = n are represented 
in Fig.l. As we can see from the figure above, the PT temperature in the f3 c A units (where A is a 




hk 

FIG. 1. The dependence of the PT temperature on the ratio of the effective range b of the interaction potential 
to the lattice constant c for different values of the n component number of the model. 

constant (see ( [2.2j )) grows with the increase of the n component number of the model. The PT temperature 
decreases and strives for the mean field theory results in the case of the increase of the effective range of 
the interaction potential b. Let us employ the introduced definitions. 



III. CALCULATION OF THE FREE ENERGY OF THE TV- VECTOR MAGNETIC MODEL. 



As shown in the previous section, the calculation of the free energy of the system can be performed in 
accodance with (2.25). The given analysis of the RR and the definition of the region of the availability of 
their approximate solutions enable us to present the free energy of the system in the vicinity of the PT 
point in the form 

F = F + Fcb + Flwf, (3.1) 
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where Fq is the free energy of the non-interacting spins 

"(27r)"/ 2 m™" 1 



Fn 



-kTNln 



r(n/2) 

Flwf is the LWF free energy, Fqr is the CR free energy 



F C R = 



Mr 

1=0 



(3.2) 



(3.3) 



where 



fi = —In 
4 



+lnU 



1 



X 



, yi-i + -r 



+ lnU 

' vt-i 



n-l 



(3.4) 



Here Ni = N's . In the case I — for /; the following relation takes place 



/ \2 




(3.5) 



The employment of the RR solutions (2.41) allows us to select in /; the explicit dependence on the number 
of the phase layer I. Having performed the summation along the layers of the CV phase space to the point 
m T of the exit of the system from the CR according to (3.3), for the free energy of the CR we obtain 

-kTN' [7q! + 7ir + 7 2 r 



F, 



CR 



■7?o 



(3.6) 



Let note that the signs " + " and " — " correspond to T > T c and T < T c respectively. The coefficients 
7oi , 71 , 72 , 7iq ar e constants and do not depend on the temperature (see Appendix 3). The analytical 
part of the CR free energy is connected with the coefficients 7qi,7i,72- It should be mentioned that the 
expressions of these coefficients coincide at the temperatures above and below the critical temperature. The 
dependence of ji on the microscopic parameters of the interaction potential and the n component number 
of the model are shown in Table III. The nonanalytical part, characterizing the temperature dependence 
of the specific heat in the vicinity of the PT point, is connected with the term jf T dL/ , where 



7l ±=7' s -^o ; 



(3.7) 



For fc R we have 



ICR 



~ d El 



2 ' 



(3.8) 



fc 



CR 



■]ny*+a(y*) 2 + 



lnf/f- 



-,**). 



(3.9) 



The values of d m are given in Appendix 4. In Table 2 the dependence of 7 10 on the component number 
n of the model for different ratios b/c is shown. The expression (3.6) describes the contribution of the 



region of the renormalization group symmetry to the free energy of the system. It allows us to obtain the 
respective contributions to the specific heat of the system at T < T c and T > T c 



C± R = kN' 



"CR 



(3.10) 



where 
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a = 2 — dv^ 

C W=2( 7l +72) ) (3.11) 
c cr = (! - a ) u cw 

U CR = dv lw- 

The curves 1 in Fig. 2 correspond to the contribution of the CR to the specific heat at T < T c and T > T c 
respectively. The negative value of the specific heat amplitude, that corresponds to the contribution of 
the CR, testifies to the nonstability of the contribution of the short-wavelength fluctuations (SWF) of the 
spin moment density. Considering the contribution of the region of the LWF of the spin moment density 
in calculating the thermodynamic characteristics of the system is a topical problem of to-day. 




■0.02 



-0.01 




-20.0 J 



-Or&h- 



-6t02 



FIG. 2. The temperature dependence of the specific heat. The comparison of the contribution of the CR and 
the region of the LWF. Notes: 1 - corresponds to the CR, 2 - limit Gaussian region (LGR) (r > 0) and inverse 
Gaussian region (IGR) (r < 0), 3 - contribution to the specific heat due to the rise of the ordering in the system. 



The regions of LGR for T > T c and IGR for T < T c correspond to the LWF of the spin moment density. 
The increase of the basic xi and the intermediate yi variables as functions of / is a typical peculiarity of the 
LGR and IGR. In this connection the contribution of these regions to the free energy can be calculated in 
the Gaussian measure density approximation (as in this region of the wave vectors the term proportional to 
the fourth power under the exponent of the distribution functions becomes much smaller than the square 
term). Since, the increase of xi is gradual, there exists the so-called transition region (TR) where it is 
necessary to keep the fourth power of the variables in the distribution function. The value of the TR 
is defined by a certain number of the CV phase space layers m following the point m T of the exit of the 
system from the CR. The value of the TR is defined by the condition 

^ (3-12) 



where a m is a constant (a m > 10). The arguments mentioned above anable us to write the contribution 
of the LWF at T > T c to the free energy of the system in the form 

FIgr = -kTN f+ GR r 3v - pp 2 B H 2 Nj 4 T- 2 » , (3.13) 

where 

/lor = Str + /*• 

It should be pointed out that the following definitions are introduced here: [ib is Bor's magneton, /tr 
corresponds to the contribution of the TR and /„ corresponds to the contribution of the wave vectors 
region with k < B s -( m r+m +i) xhe explicit analytical expressions of these coefficients was shown in 
|p0| . For 74 we have 

2m" / fn$ \2v 1 In i a\ 

74 = s fe } w (3 - 14) 

where go is a complex function of x* m , _ 1 , with rn T = m T + m +2 (sec. fl50|). 
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The contribution of the region of the LWF of the spin moment density at T < T c to the free energy of 
the system is described by the value F£ GR . The value F£ GR corresponds to the inverse Gaussian region. 
The calculation of the inverse Gaussian region contribution has its own peculiarities. In the temperature 
region T < T c the large-scale fluctuations of the spin moment density are described by the non-Gaussian 
distribution in which the coefficient near the square term becomes negative. This indicates the appearance 
of the nonzero order parameter of the system. After selecting the ordering free energy the distribution of 
the fluctuations becomes Gaussian (see Appendix 5). As a result for F£ GR we obtain 



LGR 



-kTN'\r\ 3 ^ GR , 



(3.15) 



where 



7lgr. — ^3 +73 ) 



(3.16) 



7s 



Or 



7s+7 P , 7 3 



c i/73 J 



7 3 



<S> _ „d rp/ 



The value 7-f CT> corresponds to the contribution from the ordering in the system. The coefficients j g and 
7 P are written in the form 



lg = lg c vi 



lg 



= It 



(3.17) 



3ir 2 (p n (xfj, T ) 



1 UHt 



1p = J P ct, (3.18) 
5 n ( 1 + 2f MT+ i\ n 

-% ti («i^ 2 + 3n 3 i(x)a 2 ) + 
o 

Ti 2 Q 

Ul p T + l 4/_ o-2- \ , 2- -2 

+ n (a 4 + Za^) + -n a 2 r^ +1 . 

It should be pointed out that the following definitions are introduced here 

fm T + l = /n(l + S), Um T +l = <fn ~ fnVn R l S - (3.19) 



L(x) 



arctgx 



,x 



(3.20) 



besides, c v = (j^Y is a nonuniversal value connected with the microscopic parameters of the initial 
Hamiltonian. According to (3T), (3J3), ( |3 .13 ) and ( 3.15| ) the complete expression of the free energy of the 
three dimensional n - vector magnetic model in the absence of the external field can be written as 



F = —kTN 1 [70-71H +72|r| : 



(3.21) 



where 



7o 



7oi 



sf,ln 



(27T) 



i/2 n-1 



73 ± 



r(n/2) 

= IlGR ' 



lf . 



(3.22) 
(3.23) 



The coefficient 7^ includes the contribution of the CR and the region of the LWF of the spin moment 
density at the temperatures above and below the critical. It describes the singular behaviour of the specific 
heat in the vicinity of the PT point. The dependence of the coefficients ~fo,jf on the microscopic param- 
eters of the interaction potential and the characteristics of the crystal lattice for different n component 
number of the model is shown in Tabic III. 
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IV. THERMODYNAMIC FUNCTIONS OF THE TV- VECTOR MAGNETIC MODEL IN THE 
FRAMEWORK OF THE p 4 MODEL APPROXIMATION. 



The obtained complete expression for the free energy of the n- vector magnetic model ( 3.21 ) in the p 4 
model approximation allows us to calculate ot her th ermodynamic functions in the vicinity of the phase 
transition point. As was remarked above, (see (3.10), Fig. 2), in calculating the thermodynamic functions 
of the system especially significant is considering the contribution of the LWF of the spin moment density. 
So, differentiating the expression of the free energy by the temperature, we obtain the expression for the 
entropy 



where 



S = kN [S° +S 1 T + uf\r\ 3u - 1 ] 



S°= 70 +7i, S*=2(7i +72), 



(4.1) 



The values of the coefficients S°, S 1 , uf for different n values are shown in Tabic IV 



It should be mentioned that the correct temperature behaviour of the specific heat curves is ensured 
by considering the region of the LWF of the spin moment density, i.e. the region of LGR at T > T c and 
the region of IGR at T < T c . The significant characteristic of the system is the specific heat for which we 
obtain 



where 



C = kN' C* (0) +Cf\r[ 



3i/, C<°) = 2( 7 i + 72 ), Cf = 3i/(l - a)j. 



(4.2) 



The second term in (4.2) describes the main peculiarity of the specific heat behaviour in the vicinity of the 
PT point. As we can see from (4.2), the coefficient Cf 1 includes the contributions of the CR and region of 
the LWF of the spin moment density. We can see from Fig. 2, that considering the influence of the regions 
of LWF (curves 2) provides the positivity of the specific heat and the system stability respectively. The 
dependence of the coefficients and Cf~ on the microscopic parameters of the Hamiltonian, i.e. from 
the b/c ratio for different n is exhibited in Table IV. 

In Fig. 3 the temperature dependence of the specific heat for different n is shown. As it was noticed 
above, in the case n = 3 the critical exponent a, which describes the singularity of the specific heat becomes 
negative a = -0.021 fcj, fl|. The analysis of the expression ( |4. 2[ ) and the obtained values of the specific 
heat amplitudes (see Table IV) shows that the specific heat in the case n = 3 does not diverge and receives 
a concrete finite value (see Fig. 3). 



\/5.0 -, 



10.0 : 



5.0 



0.0 



\ c 




b = c 




n= 1 










n=2 






71=3 




/ I 

- ' ' 'A 




: T 



-0.02 -0.01 QJW 0.01 



0.02 



FIG. 3. The temperature dependence of the specific heat for different n component number of the model. 

The curve of the dependence of the specific heat maximum at T = T c in the case n = 3 on the ratio 
of the effective range b of the interaction potential to the lattice constant c is shown in Fig. 4. One can 
see from this figure that the value of the specific heat maximum decreases and tends to enstant value as 
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TABLE I. The dependence of the fixed point coordinates and the eigenvalues of the transition matrix on the s 
parameter of dividing the phase space into the layers and the n component number of the model. 





s 


n 


X* 


E 1 


E2 


fn 




1.1 


1 


10.9487 


1.2077 


0.9092 


0.0181 


0.0201 




2 


11.9534 


1.2074 


0.9089 


0.0202 


0.0168 




3 


12.8801 


1.2072 


0.9086 


0.0217 


0.0144 


1.5 


1 


3.3645 


2.1521 


0.6753 


0.0998 


0.1088 




2 


3.5927 


2.1379 


0.6710 


0.1134 


0.0914 




3 


3.8088 


2.1273 


0.6676 


0.1234 


0.0787 


2.0 


1 


1.5562 


3.4761 


0.5347 


0.2153 


0.2497 




2 


1.5671 


3.3901 


0.5260 


0.2492 


0.2105 




3 


1.5848 


3.3256 


0.5185 


0.2746 


0.1814 


3.0 


1 


0.3425 


6.3985 


0.4140 


0.4640 


0.6263 




2 


0.1684 


5.9509 


0.4010 


0.5498 


0.5287 




3 


0.0154 


5.6298 


0.3880 


0.6145 


0.4538 


4.0 


1 


-0.1789 


9.6225 


0.3560 


0.7167 


1.0885 




2 


-0.4575 


8.5533 


0.3402 


0.8620 


0.9180 




3 


-0.7086 


7.8304 


0.3233 


0.9712 


0.7841 




FIG. 4. The dependence of the maximum of the specific heat on the ratio of the effective range b of the interaction 
potential to the crystal lattice constant c in the case n — 3 and T = T c . 
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TABLE II. The dependence of the temperature critical region t^ 2 on the n spin component number of the 
model. 



n 




T"2* 


1 


0.0286 


0.3449 


2 


0.0371 


0.3448 


3 


0.0424 


0.3270 



the effective range b of the interaction potential increases. This agrees with the results of the mean field 
theory. 

According to (2.3) the ratio of the critical amplitudes of the leading singular terms of the specific heat 
at T > T c and T < T c can be written in the form 

A=^L. (4.3) 

The comparison of the results obtained for the ratio of the critical amplitudes of the specific heat leading 
singular terms with results obtained by other methods is shown in Table^. 



CONCLUSIONS 



In general, it should be noted that the separate accounting of the contributions of the short- and long- 
wavelength fluctuations of the spin moment density in the expression of the free energy of the system in the 
vicinity of the PT temperature allows us to find the explicit analytical expressions for the thermodynamic 
functions as functions of the temperature. The proposed method enables as to investigate the dependence 
of the critical amplitudes of the thermodynamic functions on the microscopic parameters of the interaction 
potential and the characteristics of the crystal lattice. The results obtained for the critical exponents and 
the ratio of the critical amplitudes agree with those obtained by other methods. The negligible deviation of 
the obtained results from the experimental data and the results of the numerical calculations is connected 
with the restriction in the calculation by the p 4 model approximation. As was seen in ||, fl56| , the 
employment of the p 6 measure density for the investigation of the PT in the case n = 1 gives a more 
precise definition of the calculation results of the universal and nonuniversal characteristics of the system. 
The extension of the suggested method to the investigation of the critical behaviour of the n-vector model 
in the frame of the p 6 model approximation does not require any principal changes. To be short, in the 
present paper we do not give enough attention to calculating the scaling corrections. In the case n = 1 
such calculations were made in ||. Besides, in the common case the critical exponent of the correlation 
function is i] ^ 0. Such values of the correlation function critical exponent can be obtained by the method 
we propose if we take into account the correction on the averaging of the interaction potential in each layer 
of the CV. We are going to perform such investigations in our subsequent papers. 

APPENDIX 1 

The p"k CV are introduced by means of the functional representation for the operators of the spin density 
fluctuation p 

Pk = / PkJ (p- P 



J PwJ [p- pj (dp k ) N , p k = ■-= exp (-ikR) , 



where 



J p- p 



n 6 ~ ^) s fa ~ a*) 6 {p° ~ Po) 



is the transition operator. For the CV pk the following relations come into being 



Pk = Pk - «Pk, Pk = P^k> Pi = -F-k- 
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The Jacobian of the transition from spins to CV has the form || , J 
J[p] = Zq exp(2i7r^ p k cjk + £>M)(cKk) 

^ k 

where for Wk variables conjuagate to pk variables we have 



= g @k + i(3 k) , ^k 



W -k: ^k 



(Al.l) 



(^k) Ar =n^oii dw k cdw k s ) 



and 



Z, = 



(2tt) 2 m r ' 

rM 



(A1.2) 



For the values D[u>] the following relation is valid 



i>l ki,...k 2 



(27ri) u 2 ; ^ ,, 5 



(A1.3) 



It should be noted that here and further by the product of the vectors of the c^ki — ^k 2I type we understand 
the sum of various possible scalar products of given vectors. The coefficients U21 have the following 
expressions W2[ 



m 6m 
n n (71 + 2 J 



(A1.4) 



ii6 = 15m 6 [ 



n(n + 2)(n + 4) n 2 (n + 2) n 3 



APPENDIX 2 



For the coefficients 0,21 the following expressions are available 



«2 



n 1 
2 s o 



1 



20m 2 0$(k) 
(ri+4) JV 



C/o, 



(A2.1) 



0,4 



A S 0l 

m 4 1 



20m 2 j3$(fc) 
(n+4) AT 



T ](l-z[/„(z)-C/ 2 ). 



Here subsequent definitions are introduced 

ln + 2 
Uo = \ — ^U, 



(V), A = [B,B), 



(A2.2) 



U, 



-(«') 



■rt+l 
2 ' ' 



z = 4/ -ru 2 , 
u 4 



(A2.3) 
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where U (a, x) = D_ a _i (x) are the Weber cylinder parabolic functions. The renormalizcd u 2 i cumulants 
considering the availability of <fr(fc) for the large values of k have the form 



fceA 



u' 4 



(A2.4) 



APPENDIX 3 



For the coefficients the following expressions are valid 



A l\ , nl( A , Ai\ , _n /, , -^1 "\ I ODD. / /I I A 2 



rfe = B 7 (a 3 + ^) + 2(BiB 6 + B3B4) (a, + 
d 7 = B 5 (a 3 + ^) + 2(BiB 4 + B2B3) (a, + ^ ^ , 
d 8 - S 8 ^ 3 + ^p|) + 2 ^iB 7 + 5 2 5 6 + B3B5 + ^ (Ai + A 2 E^ 2 E^) , 

where the coefficients A t ,Bj are universal constants and are defined by the fixed point coordinates. For 
the Bj coefficients we have 

b, = ^ - ^ + \«&} , b 5 = (rir* + \ > 

B e - Vn R^-x R,-—j, B 7 - Vn — Rl ^ + — _j, 



The coefficients Ai are written in the form 

^ 1= 2 ri ~W 2 r2 ~ 4 r ' + (y*)2 ' ^3 = -2^ )' ^ = -4^V), 

where the following definition are introduced 

. (5n + 16)n n/2(n/2 + 1) ( dU n (xj) 

a = 8 2 ' n{ ~ X ) = 

r 9i 1,~ It ,3 2 1 
ri=d!-—, r 2 = -d 2 - -d iqi + - qi - -q 2 , 

U n (x*)\ dx\ ) x , ' 1 (^„(x*) V dx\ 
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APPENDIX 4 



The 7oi coefficients are expressed in the form 



7oi = s 



f* 



c 2 kd 2 E 2 c 2k diE 2 



1 - s- d 1 - s- d E? 1 - s~ d E 2 



- d ( c 2kl d 2 E 2 hd 4 E' 2 

702 = S a ( + 



703 = s 



s- d E 2 ' 1 
+ 



c\kd\E\ | c\kC 2 kd^E\E 2 c\kb§d'jE\E\ 



s~ d E 2 1 - s-^^! 1 - s- d E 1 E 2 1 - S - d E 1 E 2 ) ' 

'2 



_d / C2k2d 2 E 2 (cikc 2 ki + c-ikC\k\)d$E\E 2 c\k\d\E\ i b 2 d±E 2 



- d E 2 



1 - s- d E x E 2 



+ 



1 - s" d .Ei 



+ 



1 - 



+ 



(ci fe 6i + ci kl b )d 7 EiEl cj k d 3 E 2 c 2 lk c 2k d 6 E 2 E 2 c 2 lk b a d s E 2 E. 



1 - s- d E x El 



where 



For 7 we have 



1 - s- d E\ 1 - a- d E\E 2 1 - a- d E\El 
bo = cl k , b i = 2c 2k c 2kl , b 2 = c? 2kl + 2c 2k c 2k2 . 
f* USdi , ftt 2 d 3 
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+ 



+ 



where a definition is introduced 



1 - s- d 1 - s- d E 1 1 - s~ d Ef ' 
/„ = tE 1 t cife. 



APPENDIX 5 



The contribution of the IGR into the free energy can be written in the form 
FEgr = -kTN' S - d{m - +1 hn (2"/ 2 Q(P m j) - kTlnZ niT+1 , 



where 



Z mT+ i = / exp -i ^ d (n ' 



r + 1) 



fe<B„ 



(n,m T + l) 



4W, 



m T + l 



pit! • • • Pk 4 <5k 1 +...+k 4 ,o 



k 1 ,...k 4 



Q(P mr ) = (27T) 



-n/2 



(n,r 



-i) 



(dp) 



n/4 



<Pn(a;m T -i) 



7-7-/ ^ 1 \ / Vva T — \ \ 



(A5.1) 



(A5.2) 



After the separation of the ordering free energy in Z mr+ i the coefficient at the square term becomes 
positive. It allows us to use under the integration of Z mT+ i by the variables p k with k ^ the Gaussian 
measure density as basic. After the integration for Z mT+ i we obtain 



1 J exp {jiVN^ B Hp Q + 



7-* -Q -4 i 

B Po- jjPo f d Po- 



(A5.3) 



The value F mT+ \ corresponds to the contribution of the region of the LWF of the spin moment density 
(i.e. from the variables p\ with k — ► 0, but k ^ 0) into the free energy of the system 
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-0F, 



m T + 1 -^m T + 1 



I d (n,™ T +l)( ) | J v 



2^+1^.. v tt ; 



fc<B 

(n,mT+l) ^2 2fn,m T + l) 

°4 a 4 



(Z 4 + 3Z 2 Z 2 ) 



+ H| d (»^+i)( Q ) | 2 / 2 -^- Zx/a 

4 8 



For Zj(z = 1, 4) the next relations take place 

<,2(Mt+1) 



/x = not\- 



/3$(0) 



n 2 a 2 


■ s 2( Mt +1) 


/3$(0) 


n 3 a 3 


■ s 2(Mr+l) 


/3$(0) 


4 - 

n a4 


■ s 2(^+l) 


/3$(0) 



2 



a 2 = (a? + 6e 2 (l + e 2 )), 

/ e 3 

a 3 = a? + 6e? f 1 + -^L 



a 4 = + 6ef ( 1 + 



It should be mentioned, that 



pi 
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7T 2 [3+10f, lT + l]' 

e 2 = — |szn(7r\/2) — ttV2cos(ttV2)^ 



(A5.4) 



(A5.5) 



(A5.6) 



The variable po in (A5.3) is connected with order parameter. Its mean value is proportional to the spin 
density of the system. For the values B and G we obtain 



B=\r\ 2u B o ,B Q = p<P(0)B(°\ 
1 + S „ (0 ) 



(A5.7) 



B^ = clf n ^Bl 

E>(°) 1 " i 

Zi{ = 1 - aiiu MT+ i + a 22 



Vt+i 



2f 



G =| t r G , G = [/3«D(0)] 2 G(°), (A5.8) 



where 



«n = -- + ^" 2 a 2 , « 22 = n 3 — + . (A5.9) 

2 r Mr +i 2 I 2 3 J 

The variable po is a macroscopic value, so we can accept that 

po = VNp. (A5.10) 

It makes it possible for us to apply the saddle - point method for the calculation of Z rriT+ i in ( A5.3| ). As 
a result we find 



E>>(<p» 

where < p > is an extreme point of expression 



2tt 

, u ^^ex P {-f3F m ^ +1 ~ NE (< p >)} , (A5.ll) 



22 



E {jP> = Gi?-B? -PnHjS 



(A5.12) 



which arrises in (A5.3) with the change of variables ( A5.10). The variable po corresponds to the operator 



p Q = — = J^i cri, the mean value of which is connected with the equilibrium value of the order parameter 
a. In the case of H = we find for < p > the following solutions 



< Pi, 2 >= 




0. 



(A5.13) 



The solutions < px,2 > and < py, > correspond to extreme value of the functional Eq(p) ( A5.12 ). The 
presence of the nonzero mean spin moment at the temperatures T < T c testifies to the appearance of 
spontaneous magnetization in the system in the absence of the external field. Fig. 5 shows the temperature 
dependence of the order parameter a for different values of n component number of the model in the 
absence of the external magnetic field. 




FIG. 5. The temperature dependence of the order parameter for different values of the n component number of 
the model in the case H = (with b = c). 



In accordance with the above mentioned for Eq(< p>) we obtain 



E (<p>) = -^ = E' 



3;/ 



(A5.14) 



where 



0^2 / R (0)\2 

j Vr+i \ B i ) 



(A5.15) 



All this enables us to obtain the explicit analytical expression for the IGR free energy. Summarizing the 

-1 (fc ^ 0) and considering the contribution of the p$ variables 



expression (A5.4) by the variables k < 
in accordance with (A5.ll) - (A5.15), for F^ GR we obtain the expression (3.15) 
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TABLE III. The dependence of the coefficients of the free energy on the microscopic parameters of the interaction 
potential and the characteristics of the crystal lattice for different n component number of the model. 



b/c 


n 


7oi 


71 


72 


7r + o 


7io 


70 


7 3 + 


7 3 


).2887 


1 


0.349 


-0.500 


-0.459 


-0.538 


2.737 


1.811 


1.283 


2.726 




2 


0.727 


-0.976 


-4.427 


3.315 


0.368 


5.335 


5.795 


6.294 




3 


1.099 


-1.435 


24.497 


-25.359 


30.049 


9.488 


-22.603 


-22.530 


( = c 


1 


0.297 


-0.521 


-0.122 


-0.448 


2.276 


61.085 


1.066 


2.266 




2 


0.620 


-1.011 


-3.200 


2.722 


0.303 


192.194 


4.759 


5.170 




3 


0.938 


-1.470 


21.176 


-21.397 


25.358 


349.648 


-19.072 


-19.013 



TABLE IV. The dependence of the amplitudes of the thermodynamic functions on the microscopic parameters 
of the interaction potential and the characteristics of the crystal lattice for different n component number of the 
model. 



b/c 


n 


S° 


S 1 


«3 


u 3 


C (0) 


ct 


cr 


0.2887 


1 


1.311 


-1.917 


2.355 


-5.007 


-1.917 


1.972 


4.190 




2 


4.359 


-10.805 


11.226 


-12.196 


-10.805 


10.529 


11.436 




3 


8.053 


46.124 


-45.677 


45.530 


46.124 


-46.629 


-46.478 


b=c 


1 


60.564 


-1.286 


1.958 


-4.163 


-1.286 


1.639 


3.484 




2 


191.183 


-8.421 


9.220 


-10.017 


-8.421 


8.646 


9.393 




3 


348.178 


39.412 


-38.540 


38.422 


39.412 


-39.344 


-39.223 



TABLE V. The dependence of the ratio of the specific heat critical amplitudes at T > T c and T < T c on the n 
component number of the model. 



n 


A{This 
study) 


e — exp 
@ 


Pade— 
apr §3 


- — exp 

n 


experi- 
ment |pi[ 


1 


0.470 


0.55 


0.438 


0.519 


0.538 


2 


0.921 


0.99 


0.880 


0.888 


1.067 


3 


1.003 


1.36 


1.326 


1.152 


1.588 
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